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The use of the electrostatic classification 
method for sizing monodisperse 0.1 (im 
polystyrene latex (PSL) spheres has 
been investigated experimentally. The 
objective was to determine the feasibility 
of using electrostatic classification as a 
standard method of particle sizing in the 
development of a 0.1 (im particle diame- 
ter Standard Reference Material 
(SRM). The mean particle diameter was 
calculated from a measurement of the 
mean electrical mobility of the PSL 
spheres as an aerosol using an electro- 
static classifier. The performance of the 
classifier was investigated by measuring 
its transfer function, conducting a sensi- 
tivity analysis to verify the governing 
theoretical relationships, measuring the 
repeatability of particle sizing, and sizing 
NIST SRM 1691, 0.269 fim and NIST 
SRM 1690, 0.895 p-m particles. Investi- 
gations of the aerosol generator's perfor- 



mance focused on the effect of impuri- 
ties in the particle-suspending liquid on 
the resulting particle diameter. 

The uncertainty in particle diameter 
determined by electrical mobility mea- 
surements is found to be —3.3% to 
+ 3.0%. The major sources of uncer- 
tainty include the flow measurement, the 
slip correction, and a dependence of 
particle size on the aerosol flow rate. It 
was found that the classifier could be 
calibrated to indicate the correct size to 
within 0.1% for both SRM particle sizes 
if the defined classification length is de- 
creased by 1.9%. 

Key words: aerosol generator; atomiz- 
ers; condensation nuclei counters; 
electrical mobility; particle size; 
polystyrene latex spheres. 

Accepted: November 20, 1990 



1. Introduction 



This study assesses the acxuracy of electrostatic 
classification for measuring the diameter of 0.1 ^.m 
polystyrene latex (PSL) spheres produced as an 
aerosol by atomizing a suspension of PSL spheres 
in water. The PSL spheres used in this study were 
produced by emulsion polymerization by Dow 
Chemical Company' and the nominal size as mea- 



' Certain commercial equipment, instruments, or materials are 
identified in this paper to specify adequately the experimental 
procedure. Such identification does not imply recommendation 
or endorsement by the National Institute of Standards and 
Technology, nor does it imply that the materials or equipment 
identified are necessarily the best available for the purpose. 



sured at Dow by transmission electron microscopy 
is 0.109 \hm. This study is motivated by the need to 
develop an accurate 0.1 \im. particle size standard. 
This size standard is important for improving parti- 
cle sizing accuracy by electron microscopy, light 
scattering, and by other methods. A particle diame- 
ter of 0.1 jjtm is in the size range of combustion 
generated particles, contamination particles of con- 
cern in the semiconductor industry, air pollutant 
particulates, viruses, and various manufactured 
particulates such as carbon black and fumed silica. 
The electrostatic classifier is a widely used instru- 
ment in aerosol research for both particle sizing and 
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for generation of monodisperse aerosols over the 
size range 0.005 to 1.0 |jim. The basic physical prin- 
ciple of the classifier is that the velocity of a 
charged spherical particle in an electric field is di- 
rectly related to the diameter of the particle. A 
charged aerosol enters near the outer circumfer- 
ence of the classifier and particles with a narrow 
range in electrical mobility exit through a slit in the 
center electrode. The mobility distribution is deter- 
mined by measuring the number concentration ex- 
iting the slit as a function of the electrode voltage. 

The theory of the classifier operation and its per- 
formance have been extensively studied [1-7]. In 
regard to sizing PSL spheres with the classifier, 
Kousaka et al. [8], compared measurements of PSL 
spheres made with the electrostatic classifier to 
measurements made with a sedimentation method, 
and a balance method using a Millikan type cell. 
The three measurements were found to agree 
within a few percent for 0.2 to 1 |jim particles, but 
measurements were not conducted for particles 
smaller than 0.2 jj,m in diameter. 

The determination of the accuracy of a measure- 
ment method requires that all the physical vari- 
ables entering into the particle size equation be 
accurately known. The two key physical variables 
for the electrostatic classifier are the volumetric 
flow rate and the electrode voltage. The procedure 
used at NIST for these two calibrations is de- 
scribed in section 3.2.1. An important element in 
assessing the accuracy of an instrument is the veri- 
fication that the instrument behaves according to 
the governing equation. The verification procedure 
outlined in section 3.2.3 included comparison with 
the theory [1] and the use of two classifiers in 
tandem [9]. Further verification of the classifier 
performance is contained in section 3.2.4, where 
the measured and predicted sensitivity of the clas- 
sifier peak voltage to a change in the flow volume 
and the operating pressure are compared. An- 
other way the accuracy of the classifier was estab- 
lished was by measuring two primary calibration 
standards for particle size: 0.269 |jim (NIST SRM 
1691) and 0.895 |jim (NIST SRM 1690). The results 
of this comparison are contained in section 3.2.6. 

The method for generating the PSL sphere aero- 
sol involved atomizing a suspension of PSL spheres 
dispersed in water. The non-volatile impurities in 
PSL sphere suspension result in a residue thickness 
on the PSL sphere. Significant effort was involved 
in minimizing the droplet size produced by the at- 
omizer system, section 3.3.1, and in quantifying the 
amount of impurity in the dilution water, section 
3.3.2, and in the original, undiluted suspension, 
section 3.3.3. 



In section 4, the Discussion section, a compari- 
son is made between the results of this study and 
two other studies [10,11] that focused on the accu- 
rate measurement of the same batch of Dow 0.109 
|xm PSL spheres. 



2. Experimental Apparatus 

Figure 1. shows a schematic diagram of the in- 
strumentation used in this study. The major com- 
ponents are the atomizer, the electrostatic 
classifier, and the condensation nucleus counter. A 
PSL aerosol is produced by atomizing a suspension 
of PSL spheres in water. After conditioning, the 
particles are passed through the electrostatic classi- 
fier. By monitoring the number concentration with 
the nucleus counter versus the mobility setting of 
the classifier, the mean electrical mobility of the 
particles is determined. The mean size is then de- 
termined from the particle size dependence of the 
electrical mobility. A more detailed description of 
each instrument follows. 

2.1 Aerosol Generation 

The PSL spheres are aerosolized with an atom- 
izer, shown in figure 2, consisting of a 15 psig air jet 
impinging on the end of a liquid feed tube. The 
opposite end of the feed tube is submerged in a 
suspension of PSL spheres in water. The vacuum 
produced by the air jet draws the particle suspen- 
sion through the capillary tube and into the air jet. 
The jet atomizes the PSL particle suspension pro- 
ducing an aerosol of droplets. Some of the droplets 
produced contain PSL spheres while other droplets 
are "empty." The droplets evaporate as they flow 
through a diffusion drier and are diluted with 
clean, dry air. Droplets containing PSL spheres 
evaporate to form a PSL sphere with a slight sur- 
face residue. Droplets which do not contain a PSL 
sphere evaporate and form a small residue particle 
consisting of nonvolatile impurities present in the 
original particle suspension liquid. Thus, the result- 
ing aerosol consists of potentially dirty PSL spheres 
and small impurity particles. 

When a droplet containing a PSL sphere evapo- 
rates, any non-volatile impurities in the liquid re- 
main to form a thin layer of residue on the particle 
surface. The residue formed on the surface has a 
finite thickness and produces a systematic error in 
the measurement of particle diameter. To reduce 
the concentration of impurities in the particle 
suspension, de-ionized/filtered water was used to 
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Figure 1. Apparatus for particle sizing with the electrostatic classifier. 
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Figure 2. Atomizer with droplet impactor. 



suspend the PSL spheres. The larger droplets evap- 
orate leaving a larger impurity residue on the parti- 
cle surface. To minimize this effect, an impactor 
with a cut point of about 0.5 \xm was placed at the 
outlet of the atomizer. The effect of impurity con- 
centration on the size of 0.1 ixm PSL spheres has 
been experimentally investigated in this study. 

2.2 The Electrostatic Classifier 

The electrostatic classifier used in this project is 
a commercially available instrument (TSI, Inc., 
Model 3071). Figure 3 shows a schematic diagram 
of the instrument. The classifying region is 
bounded by a stainless steel outer cylinder with an 
inner diameter of 3.916 cm, and a coaxial, stainless 
steel center rod with an diameter of 1.874 cm. The 
center rod is connected to a variable ( to 
-11,000 V) negative dc power supply, and the 
outer cylinder is grounded. By varying the center 
rod voltage, the electric field in the annular region 
can be varied from to about 11,000 V/cm. 



Clean sheath air, after passing through a fine- 
mesh flow straightening screen at the top of the 
classifier, flows axially through the annular region 
along the center rod. A smaller, polydisperse aero- 
sol flow enters through an axisymmetric opening 
along the outer cylinder. The clean air flow forces 
the aerosol to flow downward in a thin layer on the 
outer wall of the classifying region. It is essential 
that these two streams merge smoothly without 
mixing. Near the bottom of the classifying region, a 
slit on the center rod extracts a fraction of the air 
flow consisting of near-monodisperse (single sized) 
aerosol particles. The remainder of the air flow ex- 
its through the end of the annular region as excess 
air. The length of the classifying region (44.44 cm) 
is defined as the axial distance from the aerosol 
entrance to the aerosol exit at the slit in the center 
rod. 

Before entering the classifying region, the parti- 
cles are sent through a charge neutralizer. The 
neutralization occurs through interaction with 
bi-polar gaseous ions (positive and negative ions) 
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Figure 3. Electrostatic classifier. 



produced by radioactive Kr-85. The ions impart a 
bi-polar charge distribution on the aerosol parti- 
cles. For particles with diameters around 0.1 jim, 
about 24% of the particles carry a single positive 
elementary charge and about 4% carry a double 
positive charge [12]. 

When the particles enter the classifying region, 
they are carried axially down the classifying region 
with the sheath air flow, and the particles carrying 
a positive charge move radially towards the center 



rod under the influence of the electric field. Nega- 
tively charged particles deposit on the inner sur- 
face of the outer cylinder. Within the classifying 
region, a particle rapidly reaches a steady radial 
velocity through equilibrium between the electric 
field force, and the opposing Stokes drag force. 
The radial velocity of the particle in the electric 
field is determined by the particle's electrical mo- 
bility, defined as the velocity a particle attains un- 
der the influence of a unit electric field. 
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The electrical mobility, Zp, of a singly-charged 
particle can be derived by equating the electric 
field force, Ft, with the Stokes drag force, Fi: 

Stokes Drag Force: Fi = ''L^-. , ^ 

Electric Field Force: Ft = e£ 



Electrical Mobility: 
Where 



z _v _ ^c(d;) 

^ E Sir |x Z)p 



(1) 



V = radial component of particle velocity 

E = electric field strength 

e = elementary unit of charge 

C(Dp) = slip correction 

\i = air viscosity 

Dp = particle diameter. 

As seen from eq (1), small particles have high 
electrical mobilities, and thus move with high ra- 
dial velocities toward the center rod and deposit on 
its surface. Larger particles, with lower electric mo- 
bilities, are swept further down the classifying re- 
gion before depositing on the center rod. Still 
larger particles are swept out the bottom of the 
classifier with the excess air. The monodisperse 
output of the classifier is extracted through a small 
slit on the center rod shown in figure 3. Only parti- 
cles with electric mobilities within a narrow range 
have trajectories which bring them to the entrance 
of the slit. Particles reaching the entrance of the 
slit are removed from the classifying region by the 
air flow entering the slit. In this way the classifier 
extracts a narrow size range of particles from the 
broader size range of particles entering the classify- 
ing region. 

2.2.1 The Transfer Function Knutson and 
Whitby [1] developed a theory for the classifier 
based on integrating the particle trajectory equa- 
tions. The major result of their theory is an equa- 
tion for the transfer function, H, which is defined 
as the probability an aerosol particle that enters 
the analyzer will leave via the sampling flow given 
that the particle has a mobility Zp. A brief summary 
of their analysis is presented below. 

Let r and z denote the radial and axial coordi- 
nates, respectively, within the mobility analyzer 
with z increasing in the direction of the main 
airflow. Let Ur{r;!:) and Ui{r;2) be the radial and 
axial components of the airflow velocity. Similarly, 
let Er(r,z) and Ez(r,z) be the components of the 



electric field. Neglecting particle inertia and Brow- 
nian motion, one obtains the following two first 
order differential equations for the particle path: 



dr/df = Mr -t- ZpEr, 

dz/dt = u, + ZfEz. 



(2) 
(3) 



To demonstrate the conjugate nature of the flow 
field and the electric field, Knutson and Whitby 
transformed to new coordinates, i|;, the stream func- 
tion, and <!), the electric flux function. 



r 

i|;(r,z) = [rUrdz-nizdr], 

r 

4)(r,z) = [rErdz-rEzdr]. 



(4) 



(5) 



They then demonstrate that the total differential of 
i|j-t-Zp4) is zero leading to the following significant 
result: 



v|< = —Zp4)-t- constant . 



(6) 



Quoting Knutson and Whitby [1] "The particle 
moves in such a way that the ratio of the number of 
streamlines crossed to the number of electric field 
lines crossed is always equal to the particle electric 
mobility, Zp." 

The advantage of this method of analysis is that 
the stream function is closely related to the volu- 
metric flow rate, which is an experimentally con- 
trolled variable. Representative streamlines are 
indicated in figure 4 and the corresponding flow 
variable is indicated below: 

2'ir(»]/2-v|'i) = aerosol inlet volume flow 

rate, ga 
217(11/4 - v|<2) = inlet sheath air volume flow 

rate, Qc 
217(11/4 - v|<3) = monodisperse aerosol 

volumetric flow rate, Qs 
2ir(v|i3 — t|/i) = outlet excess air volume flow 

rate, Qm. 

With the initial condition v|< = i|;in at <}> = ^i„, eq (6) 
for the particle path becomes: 

i|; = i|;i„-Zp((l)-^i„). (7) 

At <|) = <|)out, »|j has the value ^', given by: 

Vli' = lliin-ZpAcf), (8) 
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where Atf) = 4>out — <|>in. The electric field is vanish- 
ingly small in the aerosol entrance and at the exit 
slip so that (|) is a constant cj)in throughout the 
entrance and 4>out throughout the exit slit. 

The transfer function, fl, is the probability that 
the particle will leave via the sampling slit, which 
can happen only if 



l|f3<tff'<l}(4. 



(9) 



Aerosol 
Entrance 




STREAMLINES 






PARTICLE PATHS 



Figure 4. Schematic representation of mobility analyzer stream- 
lines and particle paths. 



The probability, il, is therefore equal to the fraction 
of the interval (i|fi-ZpA(|), »!i2— ZpA(|)) which is in- 
tercepted by the interval (i|;3,i|f4). The results of 
carrying out such an analysis, which is tedious but 
straightforward, are presented in figure 5. 

There are several important features of the trans- 
fer function apparent from figure 5. If the aerosol 
inlet flow and the monodisperse sampling flow are 
equal, Q^ = Q,, the transfer function has a triangular 
shape with a sharp peak corresponding to a proba- 
bility ft of 1. This is the best condition for obtaining 
accurate particle size. 
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Figure 5. The mobility analyzer transfer function. The dashed 
curve corresponds to 2s = 2a. 

For unequal flow rate, the transfer function has a 
trapezoidal shape. The origin of the top of the 
trapezoid can be understood intuitively from the 
following example. Suppose the inlet aerosol flow is 
less than the monodisperse sampling flow. Then 
there will be a range of voltages for the center rod 
for which all the inlet aerosol with mobility Zp will 
be sampled by the monodisperse outlet. This im- 
plies that the transfer function is unity for a range 
of voltages thus leading to a flat top rather than a 
triangular shaped peak. 

The actual measurements of mobility are made 
versus the voltage setting of the center rod. The 
relationship between Acj) and the voltage F can be 
obtained using eq (5) and the expression for the 
radial and axial components of the electric field: 



E,=Q,Er = VI[rln(r2lri)]. 



(10) 



Performing the integration yield the following 
result: 



A(|) = KL/ln(r2/ri). 



(11) 



The three features of figure 5 of greatest impor- 
tance to the measurement of particle size are the 
centroid of the transfer function, and the upper and 
lower widths of the transfer function. Expressing 
the results in terms of the mobility, Zp, Knutson and 
Whitby [1] obtained the following expression for the 
centroid of the mobility band. 



7 6c +6, 



fl In (rT/n) 



(12) 
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The range of the electrical mobility, AZp, of parti- 
cles exiting through the slit is given as: 



^, = ^%^Mr^n)- 



IttVL 



(13) 



The range in AZp corresponding to the upper por- 
tion of the trapezoid in figure 5, which we denote 
as AZp* is given by 



^.*=^t^^-(r^rO. 



where 



(14) 



Zp = particle electrical mobility 

AZp = electrical mobility width at base 

AZp* = electrical mobilitywidth of plateau region 

V = voltage on the center rod 

L = length from aerosol inlet to exit slit 

r2 = radius of the outer cylinder (inside surface) 

n = radius of the center rod. 

Several assumptions were made in the develop- 
ment of these equations. The flow field is assumed 
to be laminar, axisymmetric, and incompressible; 
the electric field is assumed uniform, neglecting 
field distortions at the aerosol entrance and the 
sampling exit slit; particle inertia and Brownian 
motion are neglected; and the influence of space 
and image charges are assumed negligible. 

2.2.2 Measurement of the Electrical Mobility 
Distribution By varying the voltage on the inner 
rod of the classifier, and measuring the concentra- 
tion of the aerosol exiting through the mono- 
disperse aerosol outlet, the distribution of the inlet 
aerosol's electrical mobility can be measured. The 
resolution of this measurement can be controlled, 
as seen in eqs (12) and (13), by decreasing the ratio 
of aerosol flow rate to sheath flow rate. Using eq 
(1), the electrical mobility distribution can be con- 
verted to the size distribution of the inlet aerosol. 

The classifier is instrumented with an adjustable 
voltage power supply and three mass flowmeters 
which control the sheath air, excess air, and 
monodisperse aerosol flow rates. The flowmeters 
operate by measuring the current needed to main- 
tain a constant-temperature hot-wire element in 
the air flow and are sensitive to the mass of air 
passing the sensing element. 

The calibration of a mass-sensing flowmeter can 
take the form of either an actual mass flow rate 
curve ([grams of air]/second vs meter voltage) or a 
volumetric flow rate curve ([liters of air]/second at 



T,P vs meter voltage) where T and P are the air 
temperature and pressure during calibration. As 
seen in eqs (12) and (13), measurements made with 
the classifier depend on volumetric flow rates. 
Since the classifier measures the flow rates using 
mass-sensing flowmeters, a correction must be 
made if the temperature and pressure of the air in 
the classifier differ from the temperature and pres- 
sure of the air used for the flowmeter calibration. 
If the flowmeters are calibrated using dry air, the 
correction to the calibration for dry air can be ob- 
tained from the ideal gas equation, and is as fol- 
lows: 



Qvol — Qcal ' 



cal 



I cal -t^actual 



(15) 



where 

Qvoi = volumetric flow rate 

Qcai = calibrated volumetric flow rate at TahPeai 

Teal = calibration temperature 

•Peal = calibration pressure 

Tactual = actual temperature inside classifier 

•Pactuai = actual pressure inside classifier. 

The correction to the flowmeter calibration for wet 
air is slightly different, and is described in section 
3.2.2. 

During operation of the classifier, the pressure, 
temperature and relative humidity of the air inside 
the classifier were measured, and the volumetric 
flow rate was calculated using eq (15). The temper- 
ature and relative humidity of the air inside the 
classifier were found by measuring the conditions 
of the air passing through the excess air line. The 
temperature was measured using a platinum resis- 
tance thermometer, and the relative humidity was 
measured using a chilled-mirror humidity analyzer. 
The pressure inside the classifier was about 
S.SxlCP Pa (36 cm H2O) above ambient for 
333 cmVs (20 L/min) sheath air flow. To minimize 
the gage pressure in the classifier, the excess air and 
monodisperse aerosol valves were left fully open, 
and flows were adjusted with the excess 
monodisperse aerosol valve (see fig. 3). The ele- 
vated pressure inside the classifier is required to ex- 
haust the sheath flow through the flow 
straightening header at the bottom of the classifying 
region. The pressure inside the classifier is moni- 
tored by measuring the pressure in the 
monodisperse aerosol outlet line, and applying a 
slight correction, 150 Pa (1.5 cm H2O) for a 33.3 
cmVs (2 IVmin) aerosol flow, to account for the 
pressure drop from the interior of the classifier to 
the pressure tap on the monodisperse outlet. 
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2.3 The Condensation Nucleus Counter 

Figure 6 shows a schematic diagram of the con- 
densation nucleus counter (CNC) (TSI, Inc., 
Model 3020). The instrument samples aerosol at a 
flow rate of 5 cmVs and indicates the number con- 
centration of the aerosol. The counting efficiency is 
nearly 100% for particles from about 0.02 to at 
least 0.1 [im [13]. 

The aerosol entering the counter passes through 
a chamber containing nearly saturated butyl alco- 
hol vapor. The aerosol-alcohol vapor mixture is 
then passed through a cooled condensing tube 
causing the alcohol vapor to condense onto the 



particles. The condensing alcohol causes the parti- 
cles to grow to a size easily detected with an optical 
counter at the exit of the condensing tube. In the 
optical particle counter, the particles pass through 
a focused light beam and scatter light onto a 
photodetector. In the single particle counting 
mode, used for lower particle concentrations, 
counting of individual pulses from the photodetec- 
tor provides particle concentration. In the concen- 
tration mode, used for high particle concentrations, 
the analog level of the photodetector is calibrated 
to provide particle concentration. In general, since 
the single particle counting mode does not require 
calibration, its concentration measurements are 
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Figure 6. Condensation nucleus counter. 
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considered more accurate. For sizing PSL spheres, 
the concentration was kept low enough to use the 
single particle counting mode. The PSL particle 
concentrations downstream of the classifier were 
maintained by adjusting the concentration of the 
PSL suspension used in the atomizer. 

3. Experimental Methods and Results 

Following a general description of PSL particle 
sizing with the classifier, the measurement methods 
for defining the accuracy of particle size measure- 
ments by the electrical mobility classifier are pre- 
sented. This section includes a detailed uncertainty 
analysis of the classifier performance and an anal- 
ysis of the effect of non-volatile impurities on the 
PSL sphere size as an aerosol. 

3.1 Procedure for Sizing Particles with the 
Electrostatic Classifier 

Sizing the PSL spheres with the electrostatic clas- 
sifier is a relatively fast process. A suspension of 
particles is prepared, the PSL-particle aerosol is 
generated, the classifier is used to measure the 
voltage corresponding to the mean electrical mobil- 
ity of the PSL spheres, and a straight forward data 
reduction process provides a measurement of the 
mean particle diameter. From start to finish, the 
sizing process takes about 15 min. The liquid sus- 
pension of PSL spheres was prepared by diluting a 
concentrated suspension with deionized-filtered 
water. The concentrated suspension of the 0.1 |xm 
PSL spheres consisted of about 10% by weight PSL 
spheres suspended in water. The Standard Refer- 
ence Material particles, 0.3 and 0.9 |xm, were sup- 
plied in a suspension of 0.5% by weight PSL spheres 
in water. The nominal dilutions and particle con- 
centrations of the PSL suspensions used in the 
atomizer are as follows: 



Drops of 

Particle concentrated dilution concentration 

diameter PSL suspension volume #/ml 

0.1 M.m 3 of 10% by weight 250 ml 6x10"" 

0.3 M.m 3 of 0.5% by weight 25 ml 2x10' 

0.9 n.m 10 of 0.5% by weight 25 ml 2x10" 



While the 0.3 |xm and 0.9 |xm particle concentra- 
tions in the liquid suspension were lower than the 
0.1 |xm particle concentration, the monodisperse 
aerosol concentrations were similar. The lower 



liquid concentrations of the 0.3 |xm and 0.9 |xm 
particles is offset by atomizing the suspensions 
without the impactor. The atomizer produces more 
particle-carrying droplets without the impactor. 

Following a warm-up period to allow the classi- 
fier flowmeters to stabilize, the aerosol was gener- 
ated and passed through the classifier. During 
normal operation, the sheath and excess air flow 
rates were kept equal, resulting in equal polydis- 
perse and monodisperse aerosol flow rates. To 
minimize the internal pressure of the classifier, the 
excess air and monodisperse aerosol valves were 
operated in a fully open position. The sheath air 
flow rate was set by iterating with the two valves 
upstream of the sheath air inlet, shown in figure 1, 
until the pressure between the valves was about 
1.60 X 10"* Pa (160 cm water) while maintaining the 
desired flow rate. The pressure upstream of the 
sheath air inlet was maintained at 1.60 x lO"* Pa to 
match the conditions existing during flowmeter cal- 
ibration. The excess polydisperse aerosol valve and 
the excess monodisperse aerosol valve were itera- 
tively adjusted to provide the correct excess air and 
monodisperse aerosol flow rates. The flow rates 
used for sizing the 0.1 iJim particles were nominally 
333 cmVs (20 L/min) sheath air flows, and 33.3 cmVs 
(2 L/min) aerosol flows. For sizing the 0.3 |i,m par- 
ticles, the flow rates were nominally 167 cmVs (10 
L/min) sheath flows, and 16.7 cmVs (1 L/min) aero- 
sol flows. The 0.9 |xm particles were sized using 
nominally 50 cmVs (3 LAnin) sheath, and 5 cmVs 
(0.3 L/min) aerosol flow rates. Other flow rates 
were used to investigate the effect of flow rate on 
size measurements. 

Once the flow rates in the classifier were estab- 
lished, the center rod voltage was varied to find the 
peak in the mobility distribution as measured by 
the condensation nucleus counter. The concentra- 
tion was then monitored for several minutes to 
insure a constant aerosol concentration. The fluc- 
tuations in the particle concentration were consis- 
tent with a Poisson distribution of number 
concentration; that is, the coefficient of variation, 
CV, defined as the ratio of the standard deviation 
in the number concentration to the average num- 
ber concentration, was in agreement with the pre- 
dicted CV for a Poisson distribution. 



CV = 
where 



\^ 



CV = coefficient of variation for a Poisson 

distribution 
N = average aerosol number concentration. 
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The concentration recorded for a given voltage 
setting was obtained by monitoring several consec- 
utive concentration readings (one reading every 
3 s) and then estimating the average concentration. 
When the concentration fluctuations were obvi- 
ously larger than statistically predicted, the mea- 
surement was discarded and efforts were made to 
stabilize the concentration. Gradual concentration 
changes over the course of measuring the mobility 
distribution resulted in slight sizing uncertainties 
and are included in the estimate of particle diame- 
ter measurement uncertainty. 

After the aerosol concentration stabilized, the 
voltage on the center rod of the classifier was ad- 
justed symmetrically about the peak concentration 
voltage. The concentration was recorded for each 
voltage setting. A typical concentration-voltage 
curve is shown in figure 7. The quantity of primary 
interest in this study is the peak voltage which is 
the voltage corresponding to the peak in the 
concentration-voltage curve. The peak voltage is 
computed as the concentration weighted average as 
follows: 



t^avc — 



where 






(16) 



necessary to iterate with eq (12) to determine the 
diameter. A simple iteration routine is used for this 
purpose. 

The calculation of particle diameter from eqs (1) 
and (12), requires accurate values for the viscosity 
of air, p,, and the particle slip correction, C. The 
slip correction used in the diameter calculations is 
based on Allen and Raabe's [14] measurements 
for PSL spheres using an improved Millikan 
apparatus: 



C = l-f-i&i[l.l42-f-0.558exp(^^^)l (17) 



where 



C = particle slip correction 
Kn = Knudsen number 

Kn = jr- , where Dp is the particle diameter 
\ = mean free path of air. 



Pressure and temperature corrections were 
made to the mean free path (\) [15]: 



Fave = peak voltage 

Vi = measurement voltages 

Ni = concentration corresponding to K . 

Once the representative voltage of the peak is 
found, the particle diameter can be calculated 
using eqs (1) and (12). Since the particle slip cor- 
rection is dependent on particle diameter, it is 




2000 3000 4000 5000 

Center Rod Voltage, volts 



- ^"feilpJ 



where 

Xo = 
To = 
Po = 

T = 
P = 



1-f- 



110.4 > 
To 



'1-f- 



110.4; 



(18) 



0.0673 fim, for air at To, Po 

reference temperature, 296.15 K 

reference pressure, 1.01 X 10^ Pa 

(760 mm Hg) 

air temperature; Kelvin 

air pressure inside the classifier. 



The coefficient of viscosity of air was calculated as 

[15]: 



( T \'-^/ 286.15-H10.4 \ ,,„, 

^^ = M296J5J I r-HlO.4 ) (19) 



where 



Figure 7. Number concentration vs center rod voltage for 0.1 (im 
PSL splieres. 



|X23-c= 1.93245x10-"/'. 
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3.2 Verification of Correct Classifier Performance 

To ensure correct operation of the electrostatic 
classifier, calibrations were performed and perfor- 
mance tests were conducted. The voltage and flow 
meters were calibrated, the effect of humidity on 
volumetric flow rate was measured, the transfer 
function of the classifier was measured and com- 
pared to the theoretically predicted transfer func- 
tion. A sensitivity analysis was performed to verify 
the theoretical relationships describing the classi- 
fier's dependence on the operating pressure and 
the flow rate. The repeatability of size measure- 
ments was checked and the classifier measure- 
ments were tested with Standard Reference 
Material particles. 

3.2.1 Calibration of the Flowmeters and 
Voltage Meter 

Voltage Meter Calibration 

The center rod voltage meter was calibrated with 
a precision voltage meter capable of reading 
voltages up to 10,000 V. The accuracy of the cali- 
brating meter is estimated as ±0.2%. The cali- 
bration was accomplished by connecting the cali- 
brating voltage meter to the lead from the power 
supply. The center rod voltage meter was cali- 
brated from 1000 to 9,000 V. The calibration indi- 
cated that the center rod voltage meter was 
indicating voltages higher than were actually 
present by about 2% at 4,000 and 3% at 9,000 V. 
For sizing 0.1 jjum particles using 333 cmVs sheath 
air, a 2% error in voltage at the nominal voltage 
peak of 3,800 V corresponds to a 1% error in parti- 
cle diameter. For sizing 0.3 p-m particles using 
sheath air at about 167 cmVs, a 3% error in voltage 
at a nominal voltage peak of 8,000 V corresponds 
to a 2% error in diameter. 



Flowmeter Calibration 

Calibrations of the mass flowmeters used to con- 
trol sheath air flow, excess air flow, and 
monodisperse air flow were performed to improve 
the accuracy of the size measurement. The calibra- 
tions were performed at the NIST flow calibration 
facility using the "piston prover" apparatus, main- 
tained as the primary standard for calibration of 
gas flow meters. The apparatus consists of a vol- 
ume displacement device incorporating a mercury- 
sealed piston inside of a glass cylinder. For 
calibration of a flow meter, dry gas is passed 
through the meter and into the calibration cylinder. 



The piston is displaced through an accurately de- 
fined volume in an accurately measured time. A 
bypass valve allows re-routing of the gas stream so 
the piston may be returned to its original configu- 
ration between each calibration run. Measure- 
ments of temperature and pressure are recorded so 
that the mass flow rate can be determined. 

In order to eliminate changes in the flowfield ex- 
perienced by the flowmeters, the calibrations were 
conducted without removing the flowmeters from 
the classifier. The configuration of the classifier al- 
lowed simultaneous calibration of either the sheath 
air meter and the excess air meter, or the sheath 
air meter and the monodisperse aerosol meter. To 
calibrate the sheath air meter and the excess air 
meter, the monodisperse aerosol outlet valve was 
closed and the polydisperse aerosol inlet was 
plugged. To calibrate the sheath air meter and the 
monodisperse aerosol meter, the polydisperse 
aerosol inlet was left plugged, the excess air valve 
was closed, and the monodisperse aerosol valve 
was left fully open. 

The calibrations were performed with the classi- 
fier valves in their normal configuration (excess air 
and monodisperse aerosol valves fully open). The 
valve on the sheath air inlet was adjusted to 
provide an upstream air pressure of about 
1.60x10'' Pa, and this pressure was maintained 
during normal operation of the classifier. 

During calibration, the flow rate was approxi- 
mately selected using the manufacturer's original 
calibration. The meter readings were recorded, and 
the flow rate was measured using the "piston- 
prover" calibration apparatus described above. The 
flow rates chosen for calibration of the flowmeters 
were nominally 333, 167, and 50 cmVs for the 
sheath air and excess air flow rates. These flow 
rates were chosen to maximize the flow accuracy 
for sizing 0.1 jjim, 0.3 \im. SRM, and 0.9 [xm SRM 
particles, respectively. For the monodisperse aero- 
sol flow meter, the calibration flow rates ranged 
from 33.3 to 4.2 cmVs. The flow meters were cali- 
brated at additional flow rates in the vicinity of the 
nominal values listed above. 

Each calibration point was repeated five times 
on two consecutive days, and a partial calibration 
was conducted 1 week later to check for meter 
drift. The accuracy which is normally quoted by the 
NIST calibration facility is on the order of 
±0.25%, with 99% confidence. As will be dis- 
cussed later, the estimate of uncertainty in the 
flowmeters used during operation of the classifier 
is conservatively estimated to be ±1% due to 
additional uncertainties in the meter setting and 
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the temperature and pressure which are used to 
convert the mass flow rate to the volumetric flow 
rate [eq (15)]. The calibration conducted a week 
after the initial calibration did not indicate a signif- 
icant drift for the higher flow rate calibrations 
(maximum shifts for sheath and excess air: 0.05% 
at 33.3 cmVs 0.14% at 167 cmVs, and .01% at 333 
cmVs). Drift associated with the monodisperse 
aerosol meter using lower flow rates was slightly 
higher, with the maximum shift between the three 
caUbration days of about 0.5% for flow rate set- 
tings of nominally 33.3, 16.7, and 4.2 cm^/s. 

The manufacturer's calibration for the sheath air 
meter indicated lower flow rates by about 5% at 
nominally 333 cmVs and 3% at nominally 167 cm7s 
compared to the NIST flow rate calibration. The 
manufacturer's caUbration for the excess air meter 
was found to be 8% lower at nominally 333 cmVs 
and 7% lower at nominally 167 cm^s. For sizing 0.1 
[im particles, an error in the sheath air of 5% at 
333 cm^s corresponds to a diameter error of about 
3%. It should be noted that although the electro- 
static classifier was not used until initiation of this 
project in 1988, the calibration is dated 7/83. Also, 
the larger uncertainties seen in the excess air meter 
may be due to the uncertainty in the meter setting 
caused by a significant amount of rapid fluctuation 
in the meter reading. 

3.2.2 Effect of Humidity on the Volumetric 
Flow Rate Since the molecular weight of a water 
molecule is less than the molecular weight of air, 
for a given mass flow rate, the equivalent volumet- 
ric flow rate of wet air should be higher than the 
volumetric flow rate of dry air. Water vapor, pro- 
duced by the atomizer, leads to high humidities of 
the air exiting the atomizer. While the drying tube 
and dilution of the atomizer aerosol with dry air 
reduces the humidity of aerosol entering the classi- 
fier, the resulting air humidity is still higher than 
the humidity of the air used during flowmeter cali- 
bration. The air used for sheath air is sent through 
a diffusion dryer providing relative humidities on 
the order of 5%. The relative humidity of the aero- 
sol at the classifier inlet can be high if low dilution 
air is used. The flow rate from the atomizer with- 
out dilution is 83.3 cm^/s, and a typical dilution air 
flow is about 80 cm^s. When the atomizer was used 
with the impactor, the relative humidity of the 
aerosol at the inlet to the classifier was measured 
to be around 25% and about 7% at the excess air 
outlet. 

For wet air, the volumetric flow rate correction 
made to the flowmeter calibration is slightly differ- 
ent from the correction made for dry air. Assuming 



flowmeter calibrations are conducted with dry air, 
the volumetric flow rate correction for wet air, 
derived based on ideal gas considerations, is as fol- 
lows [compare to eq (15)]: 



iSvoi =e, 



cal 



[ actual 



cal 



leal 



Pair+PHjO 



Mair 



(20) 



<2voi = volumetric flow rate of wet air 

Qcai = calibrated volumetric flow rate at Tcai/'cai 

Pair = partial pressure of air 

PhtP = partial pressure of water vapor 

Mh20 = molecular weight of water 

Mair = molecular weight of air. 

The effect of relative humidity on volumetric flow 
rate predicted by eq (20) is summarized below: 
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Typical relative humidities of the excess air mea- 
sured when sizing 0.1 nm particles were 5-15%. 
The relative humidities existing when sizing the 0.3 
and 0.9 ^m SRM particles were higher (when the 
atomizer is used without the impactor, more water 
vapor is produced). The humidities were not mea- 
sured in these cases; however, an upperbound 
humidity of 25% is estimated based on 100% 
humidity of the inlet aerosol and a factor of 10 
dilution by the dry sheath air. Particle diameter 
measurements made without correcting the flow 
rate for relative humidity will result in an increase 
in the measured diameter by a magnitude approxi- 
mately half the flow volume ratios shown above. 

The effect of humidity on the volumetric flow rate 
was experimentally investigated using a gas-test me- 
ter. Maintaining a given voltage on the mass 
flowmeter, the volumetric flow rates were measured 
with different air humidities. It was found that 
changing the humidity from 5% to 60% for fixed 
mass flow rate increased the volumetric flow rate by 
less than 0.5%, which was at the resolution limit of 
the flow measurement. This finding is consistent 
with eq (20), but the measurement resolution is 
inadequate to provide a quantitative test of the 
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equation. In any event for the 0.1 (xm PSL spheres 
with a humidity of 5-15%, the predicted humidity 
correction to the volumetric flow rate is less than 
0.2%. 

3.2.3 Testing the Transfer Function To deter- 
mine whether the classifier is operating correctly, 
its performance can be judged by comparing the 
theoretical and experimental output of the classi- 
fier when classifying a monodisperse aerosol. 
Figure 8a shows a plot of concentration vs center 
rod voltage for 0.269 |i,m PSL under the condition 
of equal aerosol flow rates. Also shown on the 
curve is the theoretical voltage-concentration curve 
plotted about the peak concentration voltage. 



which was obtained from figure 5 and eq (11) with 
Qc = Qm=167 cmVs and Ga = Gs = 17 cmVs. While 
Figure 8a indicates approximately correct behavior, 
the slight uncertainty in the peak concentration 
causes uncertainty in the placement of the theoret- 
ical transfer function. This uncertainty is the result 
of the rounding effect at the peak caused by 
slightly unequal aerosol flow rates. A better com- 
parison is obtained if the aerosol flow rates are not 
equal. The flow rate of the aerosol entering the 
classifier was reduced by a factor of two, Qa = 8.5 
cmVs, while the sheath flow was increased by 8.5 
cmVs so gc=175 cmVs (fig. 8b). The data in figure 
8b allows definite placement of the theoretical 
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Figure 8a. Comparison of experimental and theoretical response of the 
classifier for 0.269 jtm particles using equal aerosol flow rates. 
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Figure 8b. Comparison of experimental and theoretical response of the 
classifier for 0.269 ji-m particles using different aerosol flow rates. 
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transfer function and indicates correct classifier 
output. Slight differences between the theoretical 
transfer function and the experimental transfer 
function are due in part to the fact that the PSL is 
not perfectly monodisperse. 

A second method to check for correct perfor- 
mance of the classifier was developed by Rader and 
McMurry [9] and involves the use of two classifiers 
in series. Such a configuration is called a TDMA 
(Tandem Differential Mobility Analyzer). In this 
method, a polydisperse aerosol is sent through the 
first classifier to produce a test aerosol for the sec- 
ond classifier. The voltage on the first classifier is 
held constant while the voltage on the second classi- 
fier is varied to trace the distribution of the test 
aerosol. The concentration-voltage data of the sec- 
ond classifier is then compared to the TDMA 
theory using a computer algorithm which estimates 
the ratio of sheath to aerosol flow rates by fitting 
the theoretical relationships to the voltage-concen- 
tration data. Agreement between the actual flow 
ratio and the fitted flow ratio is an indication that 
the classifiers are operating correctly. This method 
was used to test the performance of the NIST clas- 
sifier using a second classifier of the same type and 
model to complete the TDMA system. The second 
classifier was provided by the University of Minne- 
sota Particle Technology Lab, The results of the 
TDMA test indicated the classifier was operating 
correctly. (For a sheath to aerosol flow ratio of 10.0, 
the algorithm indicated a ratio of 9.8 with the NIST 
classifier used as the second classifier in the TDMA 
system, iand a ratio of 10.0 with the NIST classifier 
used as the first classifier in the TDMA system.) 

A third test of the classifier's performance is to 
compare the experimental peak concentration at 
the output of the second classifier (Nzout) to the 
concentration at the input to the second classifier 
(Nia). From the triangular shape of the inlet mobil- 
ity distribution function (see fig. 5) and from a sim- 
ilar triangular shape for the sampling efficiency of 
the second classifier, Kousaka et al. [7] derived the 
following relationship between iV2out and Mn for the 
case where the voltage of the second DMA is set 
equal to the first: 



N201U = (rWin. 

(theoretical) ^^' 



(21) 
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The following experimental results indicate again 
that the classifier performs as predicted: 



The 7% to 8% difference between the actual and 
theoretical downstream concentrations is due to 
particle wall losses within the second DMA and to 
slight differences between the actual and theoreti- 
cal transfer functions. 

3.2.4 Sensitivity Analysis To investigate the 
equations governing the size measurement of the 
PSL spheres, the operating conditions of the classi- 
fier were varied slightly and then the PSL spheres 
were sized. If the equations governing particle size 
measurement are correct, the measured particle 
diameter should remain the same regardless of 
which operating conditions are used. 

The experimental method was straightforward. 
The variables which lend themselves to variation 
are flow rate and operating pressure. The flow rate 
affects the relationship between the center rod 
voltage and the particle electrical mobility 
[eq (12)], and the pressure effects the volumetric 
flow rate [eq (15)] and the particle slip correction 
[eq (17)] through its effect on the mean-free path 
of air [eq (18)]. The measurement consisted of first 
sizing 0.1 |xm PSL using nominally 333 cmVs sheath 
flow and excess flow rates, and a normal operating 
pressure of approximately 3.5 x 10^ Pa above ambi- 
ent. The particle diameter measured with these 
operating conditions was compared to the diameter 
measured when the flow rate or pressure were 
changed. Equivalently, the governing equations can 
be used to predict the change in the peak voltage 
which should result when a different flow rate or 
pressure is used for the measurement. The pre- 
dicted peak voltage can be compared to the experi- 
mentally measured peak voltage. 

Operating Pressure Variation 

The measurement using different operating 
pressures was done by restricting the excess air 
valve so that the pressure inside the classifier in- 
creased from the normal operating pressure of 
3.5x10' Pa above ambient to about 1.27x10* Pa 
above ambient. The increase in pressure results in 
a decrease in the volumetric flow rates [eq (15)], 
and a decrease in the particle slip correction factor 
[eq (17)]. The governing equations [eqs (1) and (12) 
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together with the expression for the sHp correction] 
predict that the increase in pressure should result 
in a decrease in the peak voltage of about 3%, 
which was within 0.2% of the measured decrease 
in the peak voltage, 3780 to 3685 V. The particle 
size measured with a pressure of 1.27 x 10* Pa was 
within 0.1% of the size measured using the normal 
pressure of 3.5 x 10^ Pa. The agreement between 
the two measurements of particle diameter and the 
agreement between the predicted and measured 
change in the peak voltage indicate that the pres- 
sure variable is incorporated correctly into the 
governing equations of the particle measurement. 

Sheath Flow Variation 

The sheath flow rate was decreased from 333 to 
300 cmVs while keeping the excess air flow rate at 
333 cmVs. To maintain a flow balance, the poly- 
disperse aerosol flow rate was operated at 66 cmVs, 
while the monodisperse aerosol flow rate was oper- 
ated at 33 cmVs. In this case, the governing equa- 
tions predict a decrease in the peak voltage of about 
6% which was within 0.6% of the measured change 
in peak voltage, 4240 to 4000 V. The difference be- 
tween the particle size measured using the normal 
operating conditions and the varied-flow rate condi- 
tions was less than 0.3%. This difference is probably 
caused by a dependence of measured particle 
diameter on the aerosol flow rate, which is 
discussed below. 

Aerosol Flow Variation 

Particle sizing was also conducted using different 
aerosol flow rates. In this case, the theory predicts 
a change in the mobility width of the monodisperse 
aerosol outlet particles [eq (13)], but does not pre- 
dict a change in the mean electrical mobility or 
measured particle diameter. A high aerosol flow 
rate corresponds to a wide electrical mobility range 
of the particles sampled through the slit in the cen- 
ter rod. A low aerosol flow rate corresponds to a 
narrow electrical mobility range of the mono- 
disperse aerosol outlet particles. 

To study the effect of the aerosol flow rate on 
particle size, the sheath flow and excess flow rates 
were kept constant and equal while the two aerosol 
flow rates were varied in tandem. Figures 9 and 10 
show the effect of aerosol flow rate on the voltage- 
concentration curve for 0.269 p-m and 0.1 |xm PSL 
spheres. The sheath flow and excess flow rates were 
167 cmVs for the 0.269 |xm particles and 333 cmVs 
for the 0.1 [xm particles. Figures 9 and 10 show the 
distribution plotted with both actual concentration 
and normalized concentration. Plotting normalized 



concentration allows direct comparison of the peak 
voltage. For the 0.1 |xm distribution, the particle 
diameter increases about 1% as the aerosol flow 
rate was decreased from 33.3 to 5.0 cmVs. The 
voltage peak was determined by averaging the con- 
centration-voltage data, using eq (16), for concen- 
trations greater than 0.6 A'^max. For the 0.269 |xm 
PSL spheres, an increase in diameter of about 1% 
was also found for decreasing aerosol flow rates. 

The reason for this increase in particle size is not 
presently known. The slight sizing dependence on 
aerosol flow rate is negligible for typical applica- 
tions of the electrostatic classifier. This effect 
apparently has not been reported in the literature. 
For this work, the increase in measured diameter 
for decreasing aerosol flow rates is included as an 
uncertainty in the measured diameter. 

3.2.5 Sizing Repeatability The 0.1 |xm PSL 
size measurement was repeated eight times on one 
day and six times about a week later. The 14 mea- 
surements are shown in table 1. The sheath flow 
and excess flow rates used for these measurements 
were 340 and 330 cmVs for the first and second 
days, respectively, while aerosol flow rates were 
nominally 33 cmVs. The coefficient of variation 
(CV) of the 14 measurements is 0.2%. The size 
measurements for runs 1-3 on day 1 are thought to 
have been affected by a gradually changing inlet 
concentration. If runs 1-3 on day 1 are discarded, 
the CFof the measurements is 0.1%. 

3.2.6 Measurement of Standard Reference 
Material Particles As a test of the sizing accuracy 
of the classifier. Standard Reference Material 
particles ( NIST SRM 1691 at 0.269 ±.007 |xm, and 
NIST SRM 1690 at 0.895 ±.008 |xm) were sized. 
The resulting size measurements are shown in 
table 2. The measurements were made immediately 
following the flowmeter calibration and include the 
voltage calibration. The 0.269 |xm SRM particles 
were measured using sheath flow and excess flow 
rates of 167 cmVs and aerosol flow rates of 17 cmVs. 
The 0.895 |xm SRM particles were measured using 
sheath flow and excess flow rates of 41.7 cm'/s and 
aerosol flow rates of nominally 5 cmVs. The 0.269 
ixm SRM particles, measured for the SRM report 
using electron microscopy, were measured with the 
classifier to have a mean diameter of 0.273 |xm, 
which is 1.6% larger than the SRM reported di- 
ameter. The uncertainty in the diameter of the 
0.269 |xm particles is 2.6%. The 0.895 |xm SRM 
particles, measured for the SRM report using a 
light scattering technique, were measured with the 
classifier to have a mean diameter of 0.910 |xm, 
which is 1.7% larger than the SRM reported di- 
ameter. The uncertainty in the diameter of the 
0.895 |xm particles is ±0.9%. 



162 



Volume 96, Number 2, March-April 1991 

Journal of Research of the National Institute of Standards and Technology 



400 




Le 
Symbol 


send 
Aerosol 
Flowrale 


4 


2.5 L/mIn 
1.0 UmIn 
0.5 L/mIn 



6000 7000 8000 9000 

Center Rod Voltage, volts 



Figure 9a. Voltage vs experimental concentration for 0.269 (i.m particles for three 
different aerosol flow rates. 
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Figure 9b. Voltage vs relative concentration for 0.269 ^.m particles for three 
different aerosol flow rates. 



While both measurements are larger than the 
SRM reported diameters by a similar percentage, 
the measurement of the 0.269 (jom particles lies 
within the uncertainty quoted for the SRM mea- 
surement, but the measurement of the 0.895 pom 
SRM particles is outside of the error band quoted 
for the SRM measurement. It should be noted that 
the measurement of the 0.895 ^l.m SRM particles 



was conducted using significantly different flow 
conditions (41.7 cmVs) than those used for the 0.1 
and 0.269 jim SRM particles (333 and 167 cm^s 
sheath flows, respectively). The 0.895 pom particle 
measurement can be repeated using higher sheath 
flows by measuring multiply charged particles. This 
method is described below in section 3.2.7. 
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Figure 10a. Voltage vs experimental concentration for 0.1 jim particles for three different 
aerosol flow rates. 
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Figure 10b. Voltage vs relative concentration for 0.1 jj-m particles for three different 
aerosol flow rates. 



3.2.7 Calibration of the Electrostatic ClassiHer 
Using Standard Reference Material Particles One 

possible explanation for the difference between the 
electrical mobility results for particle size and the 
certified particle size is an error in the definition of 
the length of the classifier. Recall that the length 



dimension is used in eq (12) to calculate the parti- 
cle electrical mobility from which the particle size 
is calculated using eq (1). At present, the length is 
defined as the distance from the midpoint of the 
monodisperse aerosol exit slit to the midpoint of 
the aerosol inlet (see fig. 3). This choice of length 
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Table 1. Repeatability of 0.1 Fim particle diameter measurements 



Dayl 


(May 27) 




Day 2 


(June 1) 


Run 


Measured 




Run 


Measured 


number 


diameter 
)i.m 




number 


diameter 
Fim 


1 


0.1073 




1 


0.1070 


2 


0.1072 




2 


0.1067 


3 


0.1074 




3 


0.1069 


4 


0.1070 




4 


0.1067 


5 


0.1069 




5 


0.1068 


6 


0.1069 




6 


0.1068 


7 


0.1070 








8 


0.1068 








D 


0.1071 (JLm 




D 


= 0.1068(j,m 


ffn-l = 


0.0002 (JLm 


(0.2%) 


ffn-l 


= 0.0001 fim (0.1%) 


Combined 


Analysis 








D 


0.1069 (un 








o„-i = 


0.0002 fim 


(0.2%) 







Table 2. Summary of measurements of 0.3 and 0.9 \va\ Standard Reference 
Material particles 

Standard Reference Material 1691—0.3 (j,m particles 

Measured diameter — 0.2731 (j,m 

0.2731 

0.2736 

0.2731 

_ 0.2734 

D = 0.2733 ^l.m 

o.^_, = + 0.0002 nm (0.1%) 

Dc" = 0.269 + 0.007 nm (±2.6%) 

D -De = + 0.0043 (jLm(1.6%) 



Standard Reference Material 1690—0.9 (j,m particles 

Measured diameter — 0.9103 iim 
0.9075 
0.9132 



D 


= 




0.9103 (jLm 


ffn-l 


= 




0.003 (JLm (0.3%) 


Z)c" 


= 




0.895 ± 0.008 (jLm (±0.9%) 


D -De 


= 


+ 


0.015 (JLm (1.7%) 



' Certified diameter for NIST Standard Reference Material 1691. 
'' Certified diameter for NIST Standard Reference Material 1690. 
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is consistent with the analysis by Knutson and 
Whitby [1] assuming axisymmetric and laminar 
flow and a uniform electric field in the axial direc- 
tion. These conditions will be violated to some ex- 
tent at the aerosol entrance and exit to the 
classifying column. These effects might be incorpo- 
rated in eq (12) as a corrected length of the classi- 
fier. If the length used in the calculations is taken 
as 1.9% shorter than the presently defined length, 
(44.44 cm changed to 43.60 cm), the classifier mea- 
surements of both SRM particle sizes agree within 
0.1% with the SRM reported diameters. The 
length dimension was measured in this study to be 
44.37 cm which is in close agreement with the 44.44 
cm measurement reported by the manufacturer. 
While the length can be adjusted so the classifier 
indicates the correct size for both SRM particle 
sizes, the required change in length may be too 
large to claim that the measurement differences 
are due to an error in the length definition. 

Adjusting the length definition as suggested 
above is one method of calibrating the electrostatic 
classifier for measurement of the 0.1 (j,m particles. 
A more rigorous approach for calibrating the clas- 
sifier, which is suggested for future consideration, 
is to measure the 0.269 and 0.895 [im SRM parti- 
cles using the same sheath and aerosol flow rates 
as used for the measurement of 0.1 (j,m particles. 
The calibration technique involves measuring the 
electrical mobility of multiply charged 0.269 and 
0.895 p-m SRM particles. Since a multiply charged 
particle has a higher electrical mobility than a 
singly charged particle, a higher flow rate can be 
used in the classifier to measure the mean particle 
mobility and particle diameter. By measuring the 
multiply charged SRM particles with the same flow 
conditions as the 0.1 (j,m particles, a calibration 
factor (such as changing the length definition) can 
be included in the governing equations which 
forces the SRM particle measurements to be in 
agreement with the reported diameters. This 
method of calibration is thought to be more rigor- 
ous since all the particles are measured with the 
same flow conditions. 

3.3 Investigating the EfTect of Impurities 

As was discussed in section 2.1, impurities in the 
water used to suspend the PSL spheres produce a 



layer of residue on the surface of particles after the 
PSL-cariying atomizer droplets evaporate. This 
residue results in a systematic error in particle 
diameter measurements since it increases the 
apparent particle diameter. The impurities in the 
PSL particle suspension come from impurities ex- 
isting in the water used to dilute the concentrated 
PSL particle suspension and from the impurities in 
the liquid used in the concentrated PSL particle 
suspension. To estimate the thickness of the impu- 
rity residue on the PSL sphere, it is necessary to 
know the impurity concentration in the PSL parti- 
cle suspension and the diameter of the particle-car- 
rying droplet. Assuming all of the non-volatile 
impurity forms a uniform residue shell around the 
particle, the following relationship between the 
thickness of the residue on the particle and the im- 
purity concentration, particle diameter, and 
droplet diameter is obtained: 






where 



(22) 



/= impurity addition to diameter (|xm) 

C= volumetric concentration of impurities 

Dd — PSL-carrying droplet diameter before 
evaporation ((im) 

£>p = PSL particle diameter ((im). 

To estimate the effect of impurities on the parti- 
cle size, measurements were performed to deter- 
mine the PSL particle-carrying droplet diameter 
before evaporation, the concentration of impurities 
in the water used to dilute the PSL particle suspen- 
sion, the impurities in the diluted PSL particle sus- 
pension, and the impurity concentration effect on 
particle diameter. 

3.3.1 Characterizing the Atomizer As seen in 
eq (22), the PSL-carrying droplet diameter strongly 
influences the effect of impurities on particle di- 
ameter. The droplet distribution was determined 
by atomizing a solution containing a known con- 
centration of NaCl, and measuring the resulting 
residue particle size distribution using the classifier. 
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The droplet size distribution can then be deter- 
mined using the simple relationship between the 
impurity concentration, C, residue particle size, Dp, 
and the droplet size, Dd. 



Dj, = Dd 



■ C^""^ . 



(23) 



The atomizer can be used in two configurations. 
First, for sizing the 0.1 jjim particles, an impactor is 
used at the outlet of the atomizer to remove large 
droplets (see fig. 2). Without the impactor the 
larger droplets, capable of carrying larger PSL 



spheres, are allowed to pass through the outlet of 
the atomizer. Figure 11 shows the number distribu- 
tion of droplets produced with and without the im- 
pactor. With the impactor in place, the mode of the 
number distribution is around 0.7 jjim. Without the 
impactor, the mode of the number distribution is 
around 0.8 p,m, with significantly more large 
droplets than exist with the impactor. The effect of 
the impactor is more obvious if the droplet distri- 
bution is weighted by mass or volume as shown in 
figure 12. 
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Figure 11a. Droplet number distribution with impactor. 
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Figure lib. Droplet number distribution without impactor. 



167 



Volume 96, Number 2, March-April 1991 

Journal of Research of the National Institute of Standards and Technology 



£ 


1.2 




O 






*^ 








1.0 


. 








C 






d) 






o 


O.S 


- 


o 




r 


O 






CO 
CO 

CO 


0.6 


" 






s 






<1) 


0.4 


r 


> 


















(U 


0.2 


- 


DC 


0.0 


A iiST 





10 



Dp ( n m) 
Figure 12a. Droplet mass distribution witli impactor. 




Symbol 



Legend 

NaCI 
Concentration 



530 ppm 
108 ppm 
52 ppm 



Dp(jim) 
Figure 12b. Droplet mass distribution without impactor. 



Similar measurements of an atomizer's droplet 
distribution using this technique, reported by Niida 
et al. [16], suggest that the droplet distributions 
measured in this work are biased toward larger 
particles because of diffusional losses of small 
residue particles upstream of the classifier. The 
measured droplet distributions are only qualita- 
tively representative of the actual distributions. 
The droplet distribution of the atomizer suggests 
that sizing PSL spheres without the impactor in 
place will result in larger PSL particle-carrying 
droplets, and more significant impurity effects. 

3.3.2 Measuring the Concentration of Impuri- 
ties in the Water Used to Dilute the PSL Particle 
Suspension The volumetric impurity concentra- 
tion was measured for tap water and distilled, 
deionlzed water using three methods. The tap 
water impurity concentration measurements were 
conducted for comparative purposes. Two mea- 



surement methods involved evaporating droplets 
and measuring the resulting residue particle size. 
The third method involved gravimetric measure- 
ments of evaporation residue. 

Classifying Atomized DI Water 

In the first method used to estimate water impu- 
rities, the water was atomized without PSL or other 
additives. The atomizer was used without the im- 
pactor so that larger residue particles were formed. 
The droplets formed from the atomization were 
dried and the resulting residue particles were sized 
using the classifier. This measurement was done for 
tap water and distilled, deionized (DI) water, and 
the resulting mass distributions are shown in figures 
13a and 13b along with the residue particle distribu- 
tions produced by atomizing a known solution of 
NaCl. 
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Figure 13a. Comparison of residue particles for tap water and 0.010% NaCl. 
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Figure 13b. Comparison of residue particles for deionized water and 0.0052% NaCl. 



The calculation of impurity concentration is ac- 
complished by comparing the means of the mass 
distributions of the water residue and NaCl residue 
particles. The DI water residue particles were com- 
pared to the residue particles produced from a solu- 
tion of 0.0052% NaCl by volume. The tap water 
residue particles were compared to the residue par- 
ticles produced from solution of 0.0108% NaCl by 
volume. The impurity concentrations were calcu- 
lated from the following expression derived from eq 
(22): 

Cwater = [^f^^S^\ ' X CNaC. (24) 

\X/p water mode / 

where 

Cwatet = volumetric impurity concentration 

in water 
Cnsci = volumetric concentration of NaCl 
J5pNaamode = mode of the NaCl residue distri- 
bution 
IJp water mode = mode of the Water residue distri- 
bution. 



The result of the impurity measurement follows 
from figures 13a and 13b: 



Residue 
Atomizer particle 
solution mode ((im) 


Volumetric 
impurity 
concentration 


Uncertainty 


0.0052% NaCl 0.11 
Dl-water 0.04 

0.0108% NaCl 0.15 
Tap water 0.20 


0.0052% (52 ppm) 
0.0002% (2 ppm) 
0.0108% (108 ppm) 
0.026% (260 ppm) 


± 1 ppm" 
±70 ppm 



" Uncertainties resulting from estimation of mode diameter. 

One problem with sizing the residue particles 
from the DI water is the loss of particles down- 
stream of the atomizer before being classified. To 
reduce electrostatic losses a Kr-85 bipolar charge 
neutralizer was added at the outlet of the atomizer 
for the DI water residue particles and the 0.0052% 
NaCl residue particles. Evidence of these losses is 
apparent from the observation that the number 
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concentration of residue particles for the 0.0052% 
NaCl solution is about 30 times greater than the 
concentration of the residue particles produced 
from the DI water. This would suggest that the 
mode of the number distribution of the residue 
particles from the DI water is outside the range of 
the classifier. The mode of the mass distribution of 
the residue particles appears to be within the range 
of the classifier. 

Classifying Residue Particles From the 
Vibrating Orifice 

The second method used to determine the volu- 
metric concentration of the impurity in the water 
again involved sizing the residue particles pro- 
duced by evaporating large water droplets of 
known size. In this method the Vibrating Orifice 
Monodisperse Aerosol Generator (VOAG) (TSI, 
Inc., Model 3450) was used to produce large 
monodisperse water droplets. The vibrating orifice 
generator was operated without a filter on the liq- 
uid pump so that impurities were not removed 
from the solution being tested. The VOAG was 
used to produce 39 |xm droplets of the DI water. 
The resulting residue particles were sized with the 
classifier. Except for a secondary peak correspond- 
ing to doubly-charged particles, the residue parti- 
cles were monodisperse with a size of about 0.27 
^JLm. Using eq (22), the concentration of impurities 
in the water is calculated to be 0.3 ppm. 

This same method was used to estimate the level 
of impurities in normal tap water and lab distilled 
water. Using the classifier to size the residue parti- 
cles for normal distilled water, the impurity con- 
centration was calculated to be 5 ppm. For normal 
tap water, the residual particles were too large to 
size using the classifier. Instead the TSI Model 
3310 Aerodynamic Particle Sizer (APS) was used. 
The resulting distribution is shown in figure 14 to 
have a peak particle aerodynamic diameter of 
about 2.7 |xm. For a unit density particle, the aero- 
dynamic diameter is equal to the geometric diame- 
ter, and the volumetric impurity concentration can 
be estimated using eq (22) to be 330 ppm. 

Impurity Measurements using a Gravimetric 
Method 

A third attempt at measuring the volumetric im- 
purity concentration of the water was to evaporate 
a known mass of water and measure the resulting 
impurity mass. This method did not work for the 
distilled water or the DI water because the impu- 
rity mass was too low. The method did work for the 
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Figure 14. Tap water residue particles produced by the vibrating 
orifice and measured with the aerodynamic particle sizcr. 

tap water and indicated a mass concentration of 
about 220 ppm. 

The three methods for measuring the impurity 
concentration are qualitatively in agreement. For 
the deionized/filtered water, the methods indicated 
impurity concentrations on the order of 1 ppm. The 
impurity concentration measurements of the tap 
water were conducted to compare the different 
methods. The measurements made with the 
VOAG-APS and the gravimetric method are in 
qualitative agreement, with best agreement occur- 
ring if unit density impurity is assumed. For unit 
density impurity, the VOAG-APS indicates 330 
ppm volumetric impurity and the gravimetric 
method indicating 220 ppm volumetric impurity. As 
stated above, the tap water impurity was measured 
by sizing residue particles from the atomizer to be 
about 260 ppm. 

3.3.3 Estimating the Impurity Concentration in 
the PSL Particle Suspension In the previous sec- 
tion, the impurity concentration in the water used 
to dilute. the PSL suspension was measured. The 
second source of impurities, the concentrated PSL 
suspension, is considered below. The total non- 
volatile impurity concentration in the diluted sus- 
pension is estimated in this section and also the 
predicted effect of the impurity on the particle di- 
ameter is compared with measurements. 

Calculation of the Impurity Concentration 

Concentrations of impurities in the solids com- 
posing PSL particle suspensions have been reported 
to be from about 1 to 7% [17]. If the concentrated 
PSL particle suspension is not sufficiently diluted 
with water, these impurities will have a significant 
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effect on the resulting particle size. For the 0.1 (xm 
particles, the percentage of solids in the concen- 
trated PSL suspension was about 10%. Dilution of 
the concentrated suspension was normally about 1 
to 2500 with DI water. Assuming an impurity con- 
centration in the PSL solids to be 3% (in the mid- 
dle of the range reported by Raabe [17]), the 
resulting PSL particle suspension, including 1 ppm 
impurity in the dilution water, will have an impu- 
rity concentration of about 2.2 ppm. Using eq (22) 
with a droplet diameter of about 0.9 ^,m, the re- 
sulting residue thickness (in diameter) is about 
0.00005 ^,m or about a 0.05% addition to the di- 
ameter. The same calculation for the 0.269 \im 
SRM particles containing 0.5% solids and 50 ppm 
of biocide in the concentrated PSL particle suspen- 
sion, (assuming 3% impurities in the PSL solids, a 
dilution of 1 to 250, and a larger droplet diameter 
of 2.5 [Jim since sizing occurs without the imp actor 
on the atomizer), indicates that an increase in 
diameter of 0.0001 |xm or about 0.06% would be 
expected. 



distribution of particles produced by atomizing a 
solution with a known concentration of NaCl. In 
figure 15, the entire distribution of the 0.1 \im 
PSL aerosol is shown. The secondary peaks in the 
vicinity of the main PSL particle peak are the result 
of doubly charged singlet particles and various mul- 
tiplet particles. The distribution of impurity parti- 
cles is clearly identified. The mode of the impurity 
particle mass distribution was compared to the 
mode of residue particles produced from a 
0.0052% NaCl solution indicating the impurity con- 
centration to be about 6 ppm. This same technique 
was used to estimate the impurity concentration ex- 
isting in the 0.269 SRM particles using a typical 
dilution of the concentrated PSL suspension by 
about 1 to 300 with DI water. The resulting impu- 
rity residue particles were compared to impurity 
particles from a 0.0108% NaCl solution. The mea- 
surement indicated an impurity concentration of 
about 52 ppm. This result is significantly higher 
than the estimated concentration using Raabe's 
estimates of impurities in the PSL solids. 



Measurement of the Impurity Concentration 

Measuring the concentration of impurities in the 
actual PSL suspension is difficult. As discussed in 
section 2.1, the aerosol produced by the atomizer 
consists of both PSL spheres and impurity particles 
resulting from evaporation of droplets which do 
not contain a PSL particle. Measurement of the 
impurity concentration in the PSL suspension was 
made by classifying the entire distribution of parti- 
cles existing in the PSL particle aerosol, and com- 
paring the distribution of impurity particles to the 



Effect of Impurities on the SRM Particle 
Measurements 

A simple method of determining whether or not 
the impurities in the PSL solids are contributing to 
measurement errors is to size the PSL spheres using 
a very dilute PSL particle suspension, and compare 
the measurement to that resulting from a very con- 
centrated PSL particle suspension. If the impurities 
in the undiluted PSL suspension are causing signif- 
icant measurement errors, the diameter measured 
using a very dilute suspension should be smaller 
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Figure 15. Number concentrations vs particle diameter for 0.1 p,m PSL 
particles and associated impurity residue particles. 
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than the diameter measured using a concentrated 
suspension. This measurement was conducted for 
the 0.269 (jt,m SRM particles. The concentration of 
the PSL suspension was varied by over a factor of 
15 (from 3 drops of undiluted PSL in 500 ml of DI 
water to 10 drops of undiluted PSL in 100 ml of DI 
water), and no systematic change in the diameter 
measurement was noticed. 

An attempt was also made to size the 0.269 \xm 
SRM particles using an impactor on the atomizer 
in an effort to reduce the effect of impurities in the 
suspension. The diameter measured with the im- 
pactor was similar to the diameter measured with- 
out the impactor giving strong evidence that 
impurities are not influencing the measurement. It 
is noted, however, that atomizing the 0.269 (jt,m 
particles with the impactor leads to extremely low 
particle concentrations because most of the parti- 
cles are removed by the impactor. As discussed in 
section 3.1, if the concentration of the PSL spheres 
in the aerosol is too low, significant uncertainties 
result. Based on these two measurements, it is con- 
cluded that impurities increase the particle diame- 
ter for the 0.269 p,m SRM by less than 1%. 

Effect of Impurities on Measurements of the 
0.1 \x,m Particles 

To investigate the relationship between impurity 
concentration and particle size, the nominally 0.1 
(Jim PSL spheres were sized using aqueous NaCl 
solutions of known concentration. In this measure- 
ment, the atomizer was used with the impactor in 
place because without the impactor the concentra- 
tion of NaCl particles overshadowed the PSL dis- 
tribution. The particles were first sized using clean 
water for suspension and then sized using several 
different NaCl solutions. Figure 16 summarizes the 
results. One of the data points used in the figure is 
the result obtained when the PSL is classified with 
tap water. Here the impurity concentration of the 
tap water was 0.033% as measured by sizing the 
residue particles produced by the Vibrating Orifice 
Generator with the Aerodynamic Particle Sizer. 

This data can be used to predict the effective 
PSL particle-carrying droplet diameter. The pre- 
dicted increase in diameter by eq (22) for a droplet 
diameter of 0.9 ]x,m is seen to be in good agreement 
with the data in figure 16. This diameter can be 
used, with an estimate of the impurity concentra- 
tion in the PSL suspension, to calculate the ex- 
pected residue thickness added to a particle. 
Assuming an impurity concentration in the 0.1 \x,m 
particle suspension of about 10 ppm as calculated 
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Figure 16. Impurity coneentration vs inerease in diameter for 
0.1 jLm partieles in a NaCl solution. 

and as measured, eq (22) is used to estimate the 
impurity residue increase in diameter to be about 
0.2%. Several data points taken with NaCl concen- 
trations of about 900 ppm, indicating thicknesses of 
about 0.01 \xm, were not included in Figure 16. 
This data indicated much lower thickness than 
would be expected possibly due to clumping of the 
residue on the surface of the particle. This data 
also did not agree with droplet distribution data 
obtained using lower NaCl concentrations, possibly 
because the NaCl residue particles were forming as 
hollow clumps. 

A final method for estimating the effect of impu- 
rities involved sizing of dilute and concentrated PSL 
suspensions. For PSL concentrations of 3 drops per 
1000 ml to 3 drops per 25 ml of DI water, no notice- 
able size shift occurred. This would suggest that im- 
purities are not influencing the measured diameter 
of the 0.1 ]x,m particles. If eq (22) is used with the 
previously estimated impurity concentration of 6 
ppm, and particle-carrying droplet diameter of 0.9 
]x,m, the resulting increase in diameter can be esti- 
mated to be 0.00015 ]x.m or 0.15%. The uncertainty 
estimate is approximate and we double the value 
given above so that the overall uncertainty from im- 
purities is 0.00030 ]x.m or 0.3%. 

The use of an impactor immediately downstream 
of the atomizer further reduces the impurity effect 
by removing the larger droplets. The impactor re- 
duces the peak voltage of the mobility distribution 
of the 0.1 ]x,m particles by about 130 V as indicated 
in figure 17. This corresponds to about an 0.002 |j.m 
(2%) reduction in the particle size. For the esti- 
mated impurity concentration of 6 ppm, a droplet 
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size of 2.2 jiin is estimated using eq (22). Thus it is 
found to be very important to use an impactor to 
minimize the droplet size in addition to using high 
purity dilution water. 
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Rgure 17. Effect of impactor on particle mobility. 



Estimates of Uncertainty in the ClassiOer 
Performance in the Particle Diameter 
Measurements 



Ini the previous section, results were presented 
regarding the precision associated with repeat mea- 
surements and uncertainties associated with the 
quantities appearing in the governing equations, 
eqs (1) and (12), including flow rate, voltage, and 
slip correction. In this section, a summary is pre- 
sented of all the uncertainties and an estimate of 
the overall uncertainty is given for the electrical mo- 
bility classifier. An overall estimate of the uncer- 
tainty in measuring the 0.1 \Lm PSL spheres is 
obtained by combining the uncertainty associated 
with the impurity effect and the uncertainty associ- 
ated with the use of the classifier. 



3.4,1 Random Error The random component 
of the uncertainty associated with the measure- 
ment of the average particle size can be obtained 
from the 14 repeat measurements of the particle 
size (sec. 3.2.5). The average of these 14 measure- 
ments and the associated standard deviation, a, are 
0.1068 and 0.0002 jim, respectively. The random 
component of the uncertainty, R , is given by 



R =?„-! (0.025) (r/(n) 



1/2 



(25) 



where n is the number of repeat measurements, 14, 
and t„-i (0.025) is the Student (-value forn - 1 de- 
grees of fi-eedom and for 95% confidence level 
(fi3 (0.025) = 2.16). The value of R is 0.0001 |xm, 
which corresponds to a relative error of 0.1%. 

3.4.2 Uncertainty in the Flow Rate The flow 
rate uncertainty reported in table 3 represents a 
combination of flow meter calibration accuracy, 
precision of flow rate selection, the uncertainties in 
the pressure and temperature correction to the 
flowmeter calibration, and the effect of humidity 
on the flow calibration. The precision of flow rate 
selection is estimated to be ±0.4% reflecting the 
stability of the flow rate after it is set, and the pre- 
cision of the initial flow rate setting. This value was 
calculated for the sheath air meter at nominally 
333 cmVs by estimating the precision of flowmeter 
voltage setting to be ± .002 V for a voltage setting 
of 3.160 V. The uncertainty in voltage was con- 
verted to uncertainty in flow rate using the calibra- 
tion curve. The uncertainty in the calibration, 
normally quoted by the NIST flow calibration facil- 
ity, is ±0.25% with 99% confidence. Uncertainties 
in the flow rate produced by the temperature and 
pressure correction, given by eq (15) result from 
uncertainties in the temperature and pressure. The 
uncertainty in pressure is estimated as ± 3 mm Hg 
due to uncertainties in the barometric pressure 
reading, and uncertainties in measuring the pres- 
sure inside the classifier. The uncertainty in tem- 
perature is estimated as ±0.5 °C. The resulting 
uncertainty in volumetric flow rate due to tempera- 
ture and pressure uncertainties is estimated as 
±0.4%. The effect of humidity on the volumetric 
flow rate is estimated to be 0.2% in section 3.2.2. 
The sum of all the flow related uncertainties is 
1.25% and the sum in quadrature is 0.5%. We use 
as an overall uncertainty in the volumetric flow rate 
an intermediate value of ±1.0%. 
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Table 3. Summary of uncertainties" associated with measurement of particle diameter 







Resulting 




Uncertainty 


uncertainty 


Variable 


in variable 


in diameter 


Qc = sheath air flowrate 


1.0% 


0.6% 


j2m = excess air flowrate 


1.0% 


0.6% 


rj = outer radius 


03% 


0.26% 


n = inner radius 


0.2% 


0.16% 


L = length 


0.5% 


0.30% 


V = center rod voltage 


0.45% 


0.26% 


e = elementary unit of charge 


negligible 


0.025% 


|x = viscosity of air 


0.04% 


0.025% 


C = slip correction 


0.9% 


0.5% 


T = temperature 


0.2% 


0.01% 


P = pressure 


0.4% 


0.16% 


Worst case estimate from eqs (1) and (12) 




±2.4% 


Random error, R 




±0.1% 


Residual uncertainty associated with effect 






of aerosol flowrate on apparent size 




±0.5% 


Total uncertainty associated with classifier 




±3.0% 


Impurities related uncertainty 




+0%/-0.3% 


Total uncertainty— classifier+ residue layer 




+3.0%/ -3.3% 



" The uncertainty in particle diameter determined by electrical mobility measurements 
arises from the uncertainties in the variables used in eqs (1) and (12). 



Zp = eC(Dp)/(3TT|xDp) 



Zp = (Qc+Qn,)ln(r2/n)/(4-.rl'L) 



(1) 
(12) 



3.4.3 Uncertainty in Geometric Measure- 
ments The uncertainty in the values of the center 
rod radius, n, the outer cylinder radius, ri, and the 
classification length, L, listed in table 3 are esti- 
mates of how accurately the measurements can be 
made. For the inner radius, n, the uncertainty of 
±0.2% represents about ±0.04 mm in diameter, 
which includes the variability of the diameter over 
the length of the center rod and the difference be- 
tween the diameter indicated by the manufacturer 
(1.874 cm) and the single measurement made dur- 
ing this project (1.870 cm). The uncertainty in ri is 
estimated as 0.3% in a similar manner to n, al- 
though its larger value represents the increased dif- 
ficulty in measuring the inner diameter of the 
cylinder. The uncertainty in length is estimated as 
0.5% reflecting both the uncertainty in measuring 
the length (manufacturer's measurement was 44.44 
cm compared to 44.37 cm measured in this study), 
as well as the distortion of the electric and flow 
fields at the entrance and exit of the classifier 
column. 



3.4.4 Uncertainty in Peak Voltage The uncer- 
tainty in voltage corresponding to the peak is esti- 
mated as ±0.25% corresponding to ±10 V for the 
nominal 0.1 |xm peak voltage of 3750 V when 
333 cmVs sheath air is used. This value of uncer- 
tainty was estimated by considering the data used 
during the repeatability measurement. It repre- 
sents twice the voltage spread for the second day of 
the repeatability measurements (the repeatability 
measurements are shown in table 1). The uncer- 
tainty in voltage reading due to calibration accu- 
racy is estimated as ±0.2%. Summing both 
uncertainty levels, the overall uncertainty in the 
voltage is estimated as ±0.45%. 

3.4.5 Uncertainty in Slip Correction The un- 
certainty in the slip correction for 0.1 urn particles 
is estimated to be ±0.4% based on a recent study 
by Allen and Raabe [14]. However, there is a 2.4% 
difference between the slip correction computed by 
Allen and Raabe [14] for a Knutsen number of 1.3 
( 0.1 |xm diameter sphere at ambient pressure) and 
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their earlier computation [15] based on a reanalysis 
of Millikan's oil drop data [18,19]. Allen and 
Raabe [14] attribute this discrepancy to the differ- 
ence in the surface accxjmmodation for the solid 
PSL spheres in their study compared to the liquid 
droplet surface in Millikan's studies. Because of 
this large difference (2.4%) and because there 
have been no slip correction measurements on 0.1 
|xm spheres, an intermediate estimate of the slip 
correction uncertainty of ± 0.9% is used. 

3.4.6 Pressure, Temperature, and Viscos- 
ity The effects of the uncertainty in the tempera- 
ture and pressure measurements on the flow 
calibration and flow measurements has been 
included in the flow uncertainty. However, the tem- 
perature and pressure also affect the mean free 
path of the gas eq (18), which in turn affects the slip 
correction, and the temperature affects the viscosity 
eq (19). The uncertainty in the viscosity of air itself 
is about 0.04% [20,21]. 

3.4.7 Additional Uncertainties— Aerosol Flow 
Rate and Impurity Effect All of the uncertainties 
discussed above directly affect the quantities ap- 
pearing in the governing equations. The uncertainty 
associated with varying the aerosol flow rates in 
tandem is not accounted for by propagating the un- 
certainty through the governing equations. In fact, 
as pointed out in sec. 3.2.4, the particle size com- 
puted from eqs (1) and (12) is not affected by 
changing the aerosol flow rates as long as they are 
kept equal. In section 3.2.4 the uncertainty associ- 
ated with the aerosol flow rate is estimated to be 
±0.5%. 

The uncertainty associated with impurities in the 
water does not affect the measuring accuracy of the 
classifier itself, but it does lead to a systematic in- 
crease in the PSL particle diameter as an aerosol 
compared to the actual size of the PSL sphere with- 
out any impurity coating. The impurity uncertainty 
is estimated to be 0/-0.3% based on the effect of 
impurity concentration on the residue thickness to- 
gether with the effect of an impactor on the PSL 
particle diameter (sec. 3.3.3). 

3.4.8 Total Uncertainty in the Particle Size 
Measured by the Classifier All of the various 
sources of systematic uncertainty in regard to the 
electrostatic classifier are listed in table 3. A con- 
servative measure of the combined systematic un- 
certainty is to consider the worst-case situation in 
which each variable is offset by its uncertamty to 
produce an extreme value of the diameter. The esti- 
mate is made by calculating the diameter using the 
nominal variable values and comparing to the 



diameter calculated if all the variables are offset by 
the magnitude of their uncertainty with the signs 
chosen so that the total uncertainty is a maximum. 
The percent change in the diameter is ±2.4%. 
Adding to this the random error, R, an overall 
error of ±2.5% is obtained. 

There is one additional uncertainty that must be 
included and this is the uncertainty associated with 
the aerosol flow rate, ±0.5%. Adding this value to 
the worst case total, we arrive at our best estimate 
of the uncertamty in measuring particle size with 
the electrostatic classifier as [/= ±3.0%. 

3.4.9 Total Uncertainty in the Measurement of 
the 0.1 fim PSL Spheres To obtain the total un- 
certainty in sizing the 0.1 p-m diameter PSL 
spheres with the classifier, we must include the irp- 
purity effect. While impurities in the water and in 
the particle suspension do not affect the perfor- 
mance of the classifier, they do cause the size of 
the PSL sphere to be slightly larger as an aerosol 
compared to the size of the PSL sphere itself. In 
this case the error is only in the minus direction; 
that is, this error causes the measured size to be 
too large by up to 0.3%. Adding this error to the 
worst case estimate for U given above, we obtain a 
total uncertainty for the 0.1 jjim PSL spheres of 
+ 3.0%/ -3.3%. This corresponds to the following 
range in terms of particle diameter: 



Diameter = 0.1069 



+ 0.0032 
-0.0035 



p.m. 



4. Discussion 

One way of assessing the validity of the uncer- 
tainty estimates is to compare the classifier results 
for the 0.3 and 1.0 |xm SRMs with the certified 
values. In both cases the diameter obtained by the 
classifier method is larger than the certified value, 
by 1.6% for the 0.269 |xm SRM and by 1.7% for the 
0.895 ixm SRM. The important point is that the 
percent difference between the SRM values and 
the certified values are smaller than the percent 
uncertainty ( + 3.0/ -3.3) that we have estimated 
for the 0.1 |j,m particle diameter. 

In a recent study, Knollenberg [10] summarized 
other measurements for the same batch of PSL 
sphere and reported 0.102 p,m ±0.007 |xm (Knol- 
lenberg, light scattering) and 0.105 p,m (Yamada, 
[11], electron microscopy). There are unresolved 
issues about the accuracy of size measurements by 
electron microscopy because of the uncertainties in 
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the determination of the magnification and in 
defining the edge of the particle [22], Yamada's 
study [11] has quantified the effect of the electron 
beam exposure time on the change in the particle 
diameter. The good agreement between the classi- 
fier measurements and the electron microscopy is 
encouraging but not conclusive because of the 
undefined uncertainties in the electron microscopy 
results. 

In KnoUenberg's study, the light scattering inten- 
sity of 0.1 (im PSL sphere is compared with that of 
0.269 (im SRM spheres for wavelength large 
enough that the scattering is in the Rayleigh 
regime. In this case the primary source of error is 
the uncertainty in the SRM particle itself. The size 
reported by Knollenberg [10], 0.102 jim, is outside 
the uncertainty limits of the classifier measure- 
ment; however, the uncertainty limits for the light 
scattering measurement are broad (±0.007 (im) 
and include the 0.107 jim average size obtained by 
the classifier. 

To further reduce the uncertainty associated 
with the classifier method, it is proposed that the 
classifier be calibrated with the 0.895 (im SRM, 
which has an uncertainty of ±0.9%. Both the 
0.1 (im PSL and the 0.895 (im SRM would be mea- 
sured using the same flow conditions in the classi- 
fier to remove the large flow uncertainty. By 
analyzing the multiply charged 0.895 (im particles, 
a high flow rate can be used in the classifier thus 
minimizing the uncertainties associated with oper- 
ating the classifier at low flow. It is believed that 
the uncertainty in the determination of the average 
particle size for the 0.1 jim PSL can be reduced to 
about 1.5% by using this procedure. 
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